Incommensurate Antiferromagnetism Coexisting with Superconductivity in 

Two-Dimensional d-p Model 
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Numerical studies of the two-dimensional d-p model using the Gutzwiller ansatz have exhibited 
the incommensurate antiferromagnetic state coexisting with superconductivity in the under- and 
lightly doped regions. Our results are based on the variational Monte Carlo method for the three- 
band Hubbard model with d and p orbitals. We obtained the finite superconducting condensation 
energy for the coexistent sate at the doping rate x = 1/8, 1/12, and 1/16, up to the systems of 256 
unit cells with 768 atoms (oxygen and copper atoms). The phase diagram for the hole-doped case 
is consistent with recent results reported for layered high temperature cuprates. 
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The mechanisms of superconductivity (SC) in high- 
temperature superconductors have been extensively stud- 
ied using various two-dimensional (2D) models of elec- 
tronic interactions [l], H,[!,IH]. It is of primary importance 
to clarify the phase diagram, particularly the electronic 
state in the underdoped region adjacent to the antiferro- 
magnetic (AF) phase, termed the pseudo-gap phase. It is 
unclear whether the phase diagram for L^-zSr^CuC^ is 
intrinsic for high-T c cuprates or not, although it is often 
recognized as a typical phase diagram. It is sometimes 
declared that disorder effects play some role in the spin 
glass phase of I^-zSr^CuO^. Thus, it is fair to say that 
the phase diagram has never been clarified. 

The 2D three-band d-p model is the most fundamental 
model for high-temperature cuprates U 0, 0, H, H, EH EH • 
Although we have a solution of the gap equation within 
a weak coupling perturbation theory in the limit U — » 
0[l2l[T3j. it is, however, extremely hard to show the pos- 
sibility of superconductivity exactly for finite and large 
Coulomb repulsion. Thus we adopt the Gutzwiller ansatz 
for the wave function and examine the ground state 
within the space of variational functions. We employ the 
variational Monte Carlo method [Til. [TH [l6l, [l7| to evalu- 
ate the expectation values of several physical properties. 

The purpose of this study is to investigate the co- 
existence of superconductivity and antiferromagnetism 
for the 2D d-p model. We have found that the coex- 
istent state has indeed the lowest energy in the varia- 
tional space at the doping rate x = 0.125, 0.08333, and 
0.0625 in the low-doping region. At x = 0.125, the in- 
commensurate antiferromagnetic state has eight-lattice 
periodicity, as reported on the basis of neutron scatter- 
ing measurements (l8j . The periodicity increases as x de- 
creases; we have twelve lattice periodicity at x = 0.0833 
and sixteen-lattice periodicity at x = 0.0625. 

The Hamiltonian is the d-p model containing the on- 
site Coulomb repulsion for d electrons and is written 
as 1 19] 



dp 
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d ia and d\ a are the operators for the d electrons. Pi±±/2a 
and p| ± £ / 2 <t denote the operators for the p electrons at 
the site R i±& /2, and in a similar way, p i±y /2a and p\ ± ^ /2a 
are defined. Ud is the strength of the on-site Coulomb en- 
ergy between d electrons. The number of sites is denoted 
as N s , and the total number of atoms is N a — 3N S . The 
total number of fermions is denoted as N e . The energy 
unit is given by td P in this paper. 

The van Hove singularity in the density of states plays 
an important role in two-dimensional models. We define 
the density of states as 



D(e) 



(2) 



where £k = £k — A 1 (m is the Fermi energy) and ek is the 
band crossing the Fermi energy. We examine the hole- 
doped case within the hole picture where the lowest band 
is occupied up to the Fermi energy /z. For this purpose, 
we employ the electron- hole transformation t pp —> ~t pp , 
tdp — * —tdpi and we set e p — td > 0. The density of states 
D(e) as a function of the carrier density x is shown in Fig. 
[T]for t pp = 0.4, 0.2, 0, and -0.2. x — corresponds to the 
half-filled band. For t pp = 0, the van Hove singularity is 
at x = 0. It moves to the hole-doped side of x > for 
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FIG. 1: Density of states of the d-p model as a function of 
the carrier density x for t pp = 0.4, 0.2, 0.0 and -0.2. We set 
e P = 0, ed = — 1 and td P — 1 (energy unit). 



and tpp = 0.4. Thus we set parameters to be e p — 
and t pp = 0.4 in the main computations, and Ud 



tpp > and to the electron-doped side for t pp < 0. We 
have the van Hove singularity at x ~ 0.16 for e p — = 2 

e d = 2 
= 8 in 

this paper. This is in g ood accordance with the results of 
cluster estimations [2Ct |2 Tl |22| . The van Hove singularity 
approaches x = as the level difference e p — Cd becomes 
large. Hence, we expect that the critical temperature T c 
has a peak as a function of e p ~ td if we fix the carrier 
density x. 

We adopt the Gutzwiller ansatz for the ground-state 
wave function ip: ip = PqiPo, where i/jq is a trial one-body 
wave function and 



Pg = - (! - 9>di 



n<ui) 



(3) 



is the Gutzwiller projection operator, g is the variational 
parameter in the range of < g < 1. The wave function 
considered in this paper is a coexistent state which is 
given by the solution of the Bogoliubov-de Gennes equa- 
tion: 
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j 



H 



oil 1 



(4) 



(5) 



«A= +«*H) (E~ x <0). 



(7) 



a la denotes d l(7 , Pi+x/2a, and p l+ y/ 2 a corresponding to 
the components of M^and v x . The coexistent supercon- 
ducting state is[H, 0] 



1> = PNY[a x a{\0) 



const. P/yexp 



ij 



'v a ll a ll 



10), (8) 



where |0) is the vacuum state annihilated by di a , Vi+xj2m 
and Pi+y/2a- Since ipse satisfies a\ipsc = 0, using the 
Hausdorff formula, d>a is determined as 



(9) 



where we define the matrices U and V as U\j — u x and 
V\j = Vj. Pn fixes the electron number to be N e . The 
antiferromagnetic order parameter is contained in (Hij a ) 
and the superconducting gap function is in (Fij). 

Since the incommensurate state was shown to be stable 
in the lightly doped region, we assume the spatial varia- 
tion for the order parameters. The trial Hamiltonian is 
the Hartree-Fock Hamiltonian given as[l^, [2|| [2?], [28| 

Ht„ai =K + J2[Sndi - a{-ir+y'm,]dld la . (10) 



Corresponding to the energy levels ed and e p , variational 
parameters e p and id are incorporated in the noninteract- 
ing part K in eq. (fTU)) . We assume the spatial variations 
to be 



5n di 



E 



cosh(xi — x l " c ) ' 



= A inc Y[ tanh(xj 



(11) 



(12) 



for parameters a, Ai nc , and cc™ c . x" iC determines the 



periodicity of oscillation; we set x™ c = j/(2x v ) for the 
variational parameter x v . The energy is computed for 
several values of x v such x v = 1/4, 1/8, ■ • • . A small 
spatial charge oscillation, which is, at most, ten percent 
of the total density, is induced owing to the oscillation 
potential Sridi and rrii [28[ . Thus we assume the following 
superconducting order parameter: 



A hi+X = A x cos[Q s (x l + x/2)], 



(13) 



for a trial Hamiltonian Hij a and Fij, where (Hij a ) and 
(Fij) are 3A^ S x3iV s matrices including the terms for d, p x , 
andpy orbitals. The Bogoliubov operators are written as 

a A = ^K A a iT +«,H) (^ A >0)- (6) 



Ai,i +S = A y cos(Q s Xi), (14) 

for Qs — 2irx v . We assume the d-wave symmetry for 
the SC gap function: = — A y = A. The super- 
conducting order parameter oscillates so that the am- 
plitude has a maximum in the hole-rich region and a 
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FIG. 2: Phase diagram of stable antiferromagnetic state in 
the plane of Ad p = e p — e<j and tpp obtained for 16 x 4 lattice. 
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FIG. 3: Energy of the coexistent state as a function of the 
SC order parameter for x — 0.125 on 16 x 4 lattice. We 
assume the incommensurate antiferromagnetic order (stripe) . 
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minimum in hole-poor region. The energy expectation 
value E = (ip\H\ip) / is evaluated using a Monte 

Carlo Metropolis algorithm, which is a standard method 
in variational Monte Carlo computations. 

The condensation energy E conc i is defined as the dif- 
ference E con d = E(A — ► 0) — E(A) for the optimized 
energy. The energy of the antiferromagnetic state would 
be lowered further if we consider the incommensurate 
spin correlation in the wave function. The phase di- 
agram in Fig. [5] presents the region of the stable AF 
phase in the plane of t pp and A^ p = e p — e^. For large 
A-dp = €p — ed, we have the region of the AF state with an 
eight-lattice periodicity in accordance with the results of 
neutron-scattering measurements [lg, [29(. In the incom- 
mensurate antiferromagnetic region, we obtain a finite 



FIG. 4: Energy gain due to the SC order parameter as a 
function of the system size N a tom = 3N S . Parameters are 
e p — 0, ed — —2, tpp — 0.4, and Ud = 8. The open circles 
are for the simple d-wave pairing at the hole density x = 0.2. 
The solid symbols indicate the energy gain of the coexistent 
state: the solid circles are those at x = 0.125, solid squares 
are those at x = 0.08333 and the solid triangle is that at 
x = 0.0625. The diamond shows the SC condensation energy 
obtained on the basis of specific heat measurements on the 
optimally doped YBa2Cu306+a; at x = 0.92[30l]. 
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FIG. 5: Phase diagram of the d-p model based on the 
Gutzwiller wave function. 
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SC condensation energy, assuming a spatial oscillation, 
which is shown in Fig. [3J The variational parameters are 
g = 0.386, e d = -1.578, e p = 0, fj, = -3.09, A mc = 0.5, 
and A = 0.02. 

The main results of this study are shown in Fig. [4] 
where the size dependence of the SC condensation en- 
ergy is shown for x = 0.2, 0.125, 0.08333, and 0.0625. 
We set the parameters to be e p — e d = 2 and t pp = 0.4 in 
td P units, which is reasonable from the viewpoint of the 
density of states and in the region of eight-lattice period- 
icity at x = 1/8. We have carried out the Monte Carlo 
calculations up to 16 x 16 unit cells (768 atoms in total). 
In the overdoped region in the range of 0.18 < x < 0.28, 
we have the uniform d-wave pairing state as the ground 
state. The periodicity of spatial variation judged from 
the condensation energy increases proportionally to 1 fx 
as the doping rate x decreases. In the figure, we have 
the 12-lattice periodicity at x = 0.08333 and the 16- 
lattice periodicity at x = 0.0625. For x — 0.2, 0.125, 
and 0.08333, the results strongly suggest a finite conden- 
sation energy in the bulk limit. We believe that the size 
dependence of the SC condensation energy in the incom- 
mensurate region is rather weak because the main part 
of the superfluid density is in the hole-rich region of the 
striped structure. Thus we expect a finite condensation 
energy even at x = 0.08333 and 0.0625. The SC con- 
densation energy obtained on the basis of specific heat 
measurements agrees well with the result of variational 
Monte Carlo computations [30j . In general, the Monte 
Carlo statistical errors are much larger than those for the 
single-band Hubbard model. A large number of Monte 
Carlo steps (more than 5.0xl0 7 ) is required to obtain 
convergent expectation values for each set of parameters. 

In Fig. [5] the order parameters Aaf and Asc were 
evaluated using the formula E cond = (1/2)7V(0)A 2 where 



N(0) is the density of states. The SC condensation 
energy decreases as the doping rate x is decreased be- 
cause of the striped structure of the electronic state. 
Hence, Asc also decreases. Here, we have set N(Q) ~ 
5/tdp, since N(0) is estimated to be N(0) ~ 2 to 3 
(eV)^ 1 for the optimally doped YBa 2 Cu306+2; using 
AT(0)(fc B T G ) 2 /2[3lf. The phase diagram is consistent 
with the recently reported phase diagram for layered 
cuprates [32]|. Although the incommensurate order has 
never been observed in experiments, there is a possibility 
that observed commensurability of magnetic order may 
be brought about by the effect of nearby layers, that is, 
the cancellation of incommensurability between layers. 

We examined the phase diagram of high-temperature 
superconductors with respect to the carrier density, on 
the basis of the d-p model. We carried out variational 
Monte Carlo calculations for the 2D d-p model to in- 
vestigate the ground state for large Ud- In the lightly 
doped region we obtain the coexistent state of antifer- 
romagnetism and superconductivity at the doping rate 
x = 0.125, 0.0833 and 0.0625. As long as we employ 
the Gutzwiller ansatz, the ground state exhibits coexis- 
tence in the lightly doped region. In recent experimental 
works for layered cuprates, the possibility of the coexis- 
tent state of antiferromagnetism and superconductivity 
has been ex 

plored[3l[33i. 
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